require(scales)
Loading required package: scales

Attaching package: ‘scales’

The following object is masked from ‘package:readr’:

    col_factor

Calculate averages for clusters from resolution 1.5:

v2 <- VlnPlot(seurobj, features.plot='percent.mito', point.size.use=-1, group.by='res.1.5', y.max=0.2)
v2

#Average expression is calculated by: (mean(1expm(seurobj@data[gene, cluster])))
#average.expression <- AverageExpression(SetAllIdent(seurobj, id='res.1.5'))

BATLAS results.

batlas <- read.table('../tables/tables_paper/supplementary_tables/BATLAS/BATLAS.txt', header=T, sep='\t')
batlas$cluster <- as.character(batlas$cluster)
kable(batlas)


|cluster | brown| white|
|:-------|-----:|-----:|
|0       | 0.389| 0.611|
|1       | 0.471| 0.529|
|2       | 0.384| 0.616|
|3       | 0.600| 0.400|
|4       | 0.567| 0.433|
|5       | 0.422| 0.578|
|6       | 0.315| 0.685|
|7       | 0.607| 0.393|
|8       | 0.560| 0.440|
|9       | 0.680| 0.320|
|10      | 0.479| 0.521|
|11      | 0.545| 0.455|
|12      | 0.420| 0.580|
|13      | 0.542| 0.458|
|14      | 0.723| 0.277|
|15      | 0.443| 0.557|
|16      | 0.845| 0.155|
|17      | 0.769| 0.231|
|18      | 0.488| 0.512|
|19      | 0.446| 0.554|
|20      | 0.390| 0.610|
|21      | 0.578| 0.422|

BATLAS results log normalized.

batlas <- read.table('../tables/tables_paper/supplementary_tables/BATLAS/BATLAS_log_normalized.txt', header=T, sep='\t')
batlas$cluster <- as.character(batlas$cluster)
kable(batlas)


|cluster | brown| white|
|:-------|-----:|-----:|
|0       | 0.569| 0.431|
|1       | 0.590| 0.410|
|2       | 0.564| 0.436|
|3       | 0.634| 0.366|
|4       | 0.618| 0.382|
|5       | 0.578| 0.422|
|6       | 0.539| 0.461|
|7       | 0.635| 0.365|
|8       | 0.623| 0.377|
|9       | 0.659| 0.341|
|10      | 0.594| 0.406|
|11      | 0.618| 0.382|
|12      | 0.579| 0.421|
|13      | 0.618| 0.382|
|14      | 0.677| 0.323|
|15      | 0.579| 0.421|
|16      | 0.730| 0.270|
|17      | 0.697| 0.303|
|18      | 0.595| 0.405|
|19      | 0.583| 0.417|
|20      | 0.526| 0.474|
|21      | 0.634| 0.366|
p <- ggplot(data=batlas, aes(x=cluster, y=brown)) +
  geom_bar(stat="identity") +
  ylab('Brown percentage estimate') +
  coord_flip()
p

LS0tCnRpdGxlOiAiUiBOb3RlYm9vayIKb3V0cHV0OiBodG1sX25vdGVib29rCi0tLQoKCgpgYGB7cn0KbGlicmFyeShTZXVyYXQpCmxpYnJhcnkoa25pdHIpCmxpYnJhcnkoZ2dwbG90MikKcmVxdWlyZShzY2FsZXMpCnNldXJvYmogPC0gcmVhZFJEUygnb3V0cHV0L3NldXJhdF9vYmplY3RzLzE4MDgzMS8xMHgtMTgwODMxJykKYGBgCgpDYWxjdWxhdGUgYXZlcmFnZXMgZm9yIGNsdXN0ZXJzIGZyb20gcmVzb2x1dGlvbiAxLjU6CgpgYGB7ciwgZmlnLmhlaWdodCA9IDQsIGZpZy53aWR0aCA9IDUsIGZpZy5hbGlnbiA9ICJjZW50ZXIifQpwIDwtIFRTTkVQbG90KHNldXJvYmosIGdyb3VwLmJ5PSdyZXMuMS41JywgZG8ubGFiZWw9VCwgcHQuc2l6ZT0wLjEpCnAKYGBgCgpgYGB7ciwgZmlnLmhlaWdodCA9IDEyLCBmaWcud2lkdGggPSAxMiwgZmlnLmFsaWduID0gImNlbnRlciJ9CnYxIDwtIFZsblBsb3Qoc2V1cm9iaiwgZmVhdHVyZXMucGxvdD0ncGVyY2VudC5taXRvJywgcG9pbnQuc2l6ZS51c2U9LTEsIGdyb3VwLmJ5PSdyZXMuMS41JykKdjEKYGBgCgpgYGB7ciwgZmlnLmhlaWdodCA9IDQsIGZpZy53aWR0aCA9IDEyLCBmaWcuYWxpZ24gPSAiY2VudGVyIn0KdjIgPC0gVmxuUGxvdChzZXVyb2JqLCBmZWF0dXJlcy5wbG90PSdwZXJjZW50Lm1pdG8nLCBwb2ludC5zaXplLnVzZT0tMSwgZ3JvdXAuYnk9J3Jlcy4xLjUnLCB5Lm1heD0wLjIpCnYyCmBgYAoKCmBgYHtyfQpzYXZlX3Bsb3QocCwgZmlsZT0nZmlndXJlcy9maWd1cmVzX3BhcGVyL3N1cHBsZW1lbnRhcnlfZmlndXJlcy9iYXRsYXMvMTB4LTE4MDgzMV90c25lX3Jlcy4xLjUucGRmJywgYmFzZV93aWR0aD01LCBiYXNlX2hlaWdodD00KQpzYXZlX3Bsb3QodjEsIGZpbGU9J2ZpZ3VyZXMvZmlndXJlc19wYXBlci9zdXBwbGVtZW50YXJ5X2ZpZ3VyZXMvYmF0bGFzLzEweC0xODA4MzFfdmxucGxvdF9wZXJjZW50Lm1pdG9fY2x1c3RlcmluZy4xLjUucGRmJywgYmFzZV93aWR0aD0xMiwgYmFzZV9oZWlnaHQ9MTApCnNhdmVfcGxvdCh2MiwgZmlsZT0nZmlndXJlcy9maWd1cmVzX3BhcGVyL3N1cHBsZW1lbnRhcnlfZmlndXJlcy9iYXRsYXMvMTB4LTE4MDgzMV92bG5wbG90X3BlcmNlbnQubWl0b19jbHVzdGVyaW5nLjEuNV95LWN1dG9mZi0wLjIucGRmJywgYmFzZV93aWR0aD0xMiwgYmFzZV9oZWlnaHQ9NCkKYGBgCgpgYGB7cn0KI0F2ZXJhZ2UgZXhwcmVzc2lvbiBpcyBjYWxjdWxhdGVkIGJ5OiAobWVhbigxZXhwbShzZXVyb2JqQGRhdGFbZ2VuZSwgY2x1c3Rlcl0pKSkKI2F2ZXJhZ2UuZXhwcmVzc2lvbiA8LSBBdmVyYWdlRXhwcmVzc2lvbihTZXRBbGxJZGVudChzZXVyb2JqLCBpZD0ncmVzLjEuNScpKQpgYGAKCkJBVExBUyByZXN1bHRzLgoKYGBge3J9CmJhdGxhcyA8LSByZWFkLnRhYmxlKCd0YWJsZXMvdGFibGVzX3BhcGVyL3N1cHBsZW1lbnRhcnlfdGFibGVzL0JBVExBUy9CQVRMQVMudHh0JywgaGVhZGVyPVQsIHNlcD0nXHQnKQpiYXRsYXMkY2x1c3RlciA8LSBhcy5jaGFyYWN0ZXIoYmF0bGFzJGNsdXN0ZXIpCmthYmxlKGJhdGxhcykKYGBgCgpgYGB7ciwgZmlnLmhlaWdodCA9IDYsIGZpZy53aWR0aCA9IDEyLCBmaWcuYWxpZ24gPSAiY2VudGVyIn0KcCA8LSBnZ3Bsb3QoZGF0YT1iYXRsYXMsIGFlcyh4PWNsdXN0ZXIsIHk9YnJvd24pKSArCiAgZ2VvbV9iYXIoc3RhdD0iaWRlbnRpdHkiKSArCiAgeWxhYignQnJvd24gcGVyY2VudGFnZSBlc3RpbWF0ZScpICsKICBjb29yZF9mbGlwKCkKcApgYGAKCmBgYHtyfQpzYXZlX3Bsb3QoJ2ZpZ3VyZXMvZmlndXJlc19wYXBlci9zdXBwbGVtZW50YXJ5X2ZpZ3VyZXMvYmF0bGFzL0JBVExBU19yZXN1bHRzLnBkZicsIHAsIGJhc2VfaGVpZ2h0PTYsIGJhc2Vfd2lkdGg9MTIpCmBgYAoKQkFUTEFTIHJlc3VsdHMgbG9nIG5vcm1hbGl6ZWQuCgpgYGB7cn0KYmF0bGFzIDwtIHJlYWQudGFibGUoJ3RhYmxlcy90YWJsZXNfcGFwZXIvc3VwcGxlbWVudGFyeV90YWJsZXMvQkFUTEFTL0JBVExBU19sb2dfbm9ybWFsaXplZC50eHQnLCBoZWFkZXI9VCwgc2VwPSdcdCcpCmJhdGxhcyRjbHVzdGVyIDwtIGFzLmNoYXJhY3RlcihiYXRsYXMkY2x1c3RlcikKa2FibGUoYmF0bGFzKQpgYGAKCmBgYHtyLCBmaWcuaGVpZ2h0ID0gNiwgZmlnLndpZHRoID0gMTIsIGZpZy5hbGlnbiA9ICJjZW50ZXIifQpwIDwtIGdncGxvdChkYXRhPWJhdGxhcywgYWVzKHg9Y2x1c3RlciwgeT1icm93bikpICsKICBnZW9tX2JhcihzdGF0PSJpZGVudGl0eSIpICsKICB5bGFiKCdCcm93biBwZXJjZW50YWdlIGVzdGltYXRlJykgKwogIGNvb3JkX2ZsaXAoKQpwCmBgYAoKYGBge3J9CnNhdmVfcGxvdCgnZmlndXJlcy9maWd1cmVzX3BhcGVyL3N1cHBsZW1lbnRhcnlfZmlndXJlcy9iYXRsYXMvQkFUTEFTX3Jlc3VsdHNfbG9nLW5vcm1hbGl6ZWQucGRmJywgcCwgYmFzZV9oZWlnaHQ9NiwgYmFzZV93aWR0aD0xMikKYGBgCgoK